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SECTION  I 


INTRODUCTION 

The  NonHectrical  Languages  Simulation  Module  (NELSIM)  for  the  Air 
Force  Weapons  Laboratory  Systems  Analysis  Program  provides  the  capability 
for  large-scale  nonlinear  system  analysis  of  systems  composed  of  mechanical, 
thermal,  electromechanical,  and  electro-optical  components  in  addition  to 
electrical  and  electronic  components.  This  module  permits  engineers  not 
familiar  with  disciplines  generally  associated  with  circuit  simulation 
programs  to  solve  system  problems  without  manually  deriving  electrical 
analogs  and  differential  equations  for  use  with  the  circuit  analysis  pro- 
grams . 

A number  of  computer-aided  circuit  analysis  programs  have  been 
developed  to  automatically  calculate  the  transient  response  of  nonlinear 
circuits.  To  use  these  programs,  it  is  generally  necessary  to  describe 
the  circuit  topology  to  the  program  in  terms  of  the  electrical  network 
parameters  resistance,  inductance,  capacitance,  voltage  and  current.  The 
programs  automatically  generate  a set  of  algebraic  and  differential 
equations  from  the  network  topology;  the  algebraic  equations  are  solved 
and  the  differential  equations  are  numerically  integrated  to  determine 
the  transient  response. 

Current  and  future  weapons  systems  have  rigid  specifications  for 
operation  in  benign  and  hostile  nuclear  weapon  environments.  It  has 
become  increasingly  important  to  have  efficient  computer-aided  system 
analysis  capabilities  in  addition  to  the  existing  circuit  analysis 
capabilities.  The  system  problem  can  be  solved  using  numerical  methods 
similar  to  those  used  in  circuit  analysis  programs  because  systems, 
like  circuits,  can  be  defined  in  terms  of  algebraic  and  differential 
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equations.  Therefore,  it  is  possible  to  adapt  circuit  analysis  programs 
to  the  solution  of  system  level  problems. 

While  the  system  to  be  analyzed  can  be  mathematically  described 
in  a manner  consistent  with  available  circuit  analysis  programs,  the 
systems  may  contain  components  other  than  electrical  or  electronic 
components,  and  to  predict  the  transient  system  response,  it  is  necessary 
to  include  the  non-electrical  components  in  the  system  analysis.  In 
adopting  the  circuit  analysis  programs  to  system  level  problems,  an  auto- 
matic language  processor  was  developed  to  convert  descriptions  of  non- 
electronic system  components  to  a form  compatible  with  the  applicable 
analysis  programs  being  used.  This  procedure  permits  the  description  of 
the  entire  system  to  the  program  in  the  language  particular  to  the  differ- 
ent system  components.  The  language  processor  automatically  develops  the 
appropriate  algebraic  and  differential  equations  to  determine  the  transient 
response;  and  the  system  analyst  is  relieved  of  developing  electrical  ana- 
logs or  writing  differential  equations  compatible  with  the  system 
simulation  code  of  interest. 

The  program  translates  the  given  system  inputs  into  electrical 
analogues  and  equations  compatible  with  the  SCEPTRE  program.  To  give  the 
program  flexibility  to  interface  with  the  AFWL  systems  analysis  program 
and  provide  the  capability  to  interface  with  transient  analysis  programs 
other  than  SCEPTRE,  the  program  is  divided  into  three  main  parts  - an 
input  processor,  a translator,  and  an  output  processor  as  shown  in 
Figure  1 . 
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Figure  1.  Program  Modular  Configuration 

The  input  processor  interprets  and  stores  all  input  information. 

This  input  can  be  in  three  forms,  illustrated  in  Figure  2.  The  trans- 
lator generates  the  electrical  analogs  based  upon  the  information  stored 
within  the  program.  The  output  processor  then  generates  analogous  in  a 
format  acceptable  to  SCEPTRE.  To  alter  the  program  to  be  compatible  with  '' 

other  programs  simply  involves  alteration  of  the  output  module. 

The  remainder  of  this  report  documents  the  theory  and  formulation 
utilized  in  the  generation  of  NELSIM.  Sections  II  through  IV  document 
those  portions  of  the  program  dealing  with  analog  system  functional 
elements,  differential  equations  and  transfer  functions,  respecti vely. 

Section  V documents  the  program  configuration  and  functional  flow.  The 
versatility  of  the  program  is  illustrated  through  the  sample  problem 
package  contained  in  Section  VI. 
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Figure  2.  Input/Output  of  NELSIM 
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SECTION  II 


CONVERSION  TO  ELECTRICAL  ANALOGS 

Systems  of  equations  from  one  scientific  discipline  that  have  the 
same  form  as  those  equations  of  another  discipline  are  called  analogs. 

The  quantities  which  occupy  corresponding  positions  in  the  two  equations 
are  called  analogous  quantities.  This  section  is  concerned  with  establish- 
ing analogs  from  the  non-electrical  disciplines  to  the  electrical  disci- 
pline using  NELSIM.  Electrical  analogs  allow  solution  of  non-electrical 
problems  using  computer  programs  such  as  SCEPTRE,  NET-2  or  CIRCUS  II. 

The  analogs  from  mechanical  and  thermal  systems  to  electrical  systems 
are  summarized  in  Table  I. 

Table  I.  Analogous  Quantities  Between  Mechanical, 

Thermal  and  Electrical  Systems 


generalized 

QUANTITIES 

ELECTRICAL 

QUANTITIES 

MECHANICAL 

translational 

quantities 

MECHANICAL 

ROTATIONAL 

QUANTITIES 

THERMAL 

QUANTITIES 

1 

Flow  Variable 

q(t) 

Current 

i(t) 

Force 

f(t) 

Torque 

T{t) 

Heat  Flow 
Rate 

q(t) 

Potential 

Variable 

r(t) 

Vol tage 
e(t) 

Velocity 

v(t) 

Angular 

Velocity 

^,(t) 

Temperature 

9(t) 

POWER 

e(t)i{t) 

nt)v(t) 

T(t)^^(t) 

q(t) 

q(t)  » Gp(t) 

i(r)  = ^(t) 

Ht)  = Bv(t) 

i(t)  = B^oi^(t) 

q(t)  = 

i(t)  = C 

f(t)=M^ 

(t)  = 

q(t)  = C, 

q(t)  * K p dt 

i(t)  = i y"e  dt 

f(t)  = kJv  dt 

r(t)  = 
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The  disciplines  allowed  as  input  into  the  NELSIM  program  are  listed  below 
and  discussed  individually  in  the  following  paragraphs. 

1.  Mechanical 

2.  Thermal 

3.  Electro-mechanical 

4.  Electro-optical 

a.  Mechanical/Elettrical  Analogs 

Mechanical  systems  can  be  described  by  various  combinations  of  trans- 
lational and  torsional  elements  and  forces.  The  analysis  of  dynamic 
mechanical  systems  primarily  consists  of  determining  the  time  dependence 
of  resulting  translational  or  torsional  displacements,  velocities,  and 
accelerations  of  the  various  system  elements  under  the  influence  of  an 
externally  applied  disturbing  force  or  torque.  The  basic  description  of 
a mechanical  system  includes  the  values  of  the  system  lumped  elements  and 
their  physical  interconnection  with  eacn  other.  This  is  in  addition  to 
specifying  the  magnitude  and  time  behavior  of  the  externally  applied  dis- 
turbing force  or  torque.  The  dissipating  or  damping  parameter  is  the  dash- 
pot  or  friction  device.  For  linear  friction  or  damping  models  the  retard- 
ing force  is  proportional  to  the  velocity  and  in  a direction  opposite  to  the 
velocity  vector.  Friction  of  this  linear  nature  is  known  as  viscous 
friction,  or  viscous  damping.  Its  existence  is  well  established  experi- 
mentally for  moderate  velocities.  For  large  velocities,  the  resistance 
becomes  nonlinear  and  may  be  more  nearly  proportional  to  the  square  or  even 
the  cube  of  the  velocity. 

The  mass,  moment  of  inertia,  spring  and  elastic  shaft  act  as  energy 
reservoirs  of  the  system.  The  spring  and  elastic  shaft  are  reservoirs  of 
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potential  energy  and  the  mass  and  moment  of  inertia  are  reservoirs  of 
kinetic  energy.  The  energy  storing  characteristic  of  the  spring  is 
expressed  in  terms  of  the  force  per  unit  translational  displacement. 
Similarly,  the  elastic  shaft  has  units  of  torque  per  unit  angular  dis- 
placement. A nonlinear  spring  would  exist  for  the  condition  when  the 
spring  is  stretched  beyond  its  elastic  limits  and  the  plastic  deformation 
of  the  material  would  no  longer  produce  a spring  displacement  proportional 
to  the  applied  force. 

The  energy  storing  characteristic  of  the  mass  is  related  in  terms  of 
the  force  per  unit  acceleration  or  weight  per  unit  gravitational  accelera- 
tion. Similarly,  the  moment  of  inertia  has  units  of  torque  per  unit 
angular  acceleration.  An  example  of  non-constant  mass  situation  which 
would  result  in  nonlinear  differential  equations  arises  in  relativistic 
mechanics  where  the  effective  mass  of  a particle  increases  as  the  velocity 
of  the  particle  approaches  the  speed  of  light.  An  infinite  force  is 
required  to  propel  a particle  with  even  the  smallest  rest  mass  at  the 
speed  of  light. 

For  a mechanical  translational  system,  the  dynamic  analysis  of  the 
system  is  based  upon  Newton's  law: 

Mass  X Acceleration  = Force 

The  analysis  of  a torsional  system  is  based  upon  Newton's  law  in  the  form: 
Moment  of  Inertia  x Angular  Acceleration  = Torque 
The  elements  involved  in  translational  and  rotational  systems  are  illus- 
trated in  Figure  3 along  with  their  respective  relationships.  Table  II 
shows  the  units  involved  for  each  type  of  system.  A set  of  equations 
describing  the  dynamic  behavior  of  a mechanical  system  is  obtained  by 
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Table  II.  Mechanical  Units 


TRANSLATIONAL 


ROTATIONAL 


fit) 

v{t) 

M 

B 

K 

mks  units 

Newton 

meter/sec 

Kilogram 

(kg) 

Newton-sec/ 

meter 

Newton/ 

meter 

English 

Units 

pound  (1b) 

ft/sec 

slug 

Ib-sec/ft 

lb/ ft 

T(t) 

Wp(t) 

J 

«r 

V 

mks  units 

Newton- 

meter 

radians/sec 

Kilogram- 

meter^ 

Newton- 

meter 

sec/radian 

Newton - 

meter 

radian 

English 

units 

Ib-ft 

radians/sec 

slug-ft^ 

Ib-ft-sec/ 

radian 

lb- ft/ 
radian 
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setting  the  sum  of  the  reacting  forces  (or  torques)  equal  to  the  disturbing 
or  applied  force  (or  torque).  This  is  equivalent  to  setting  all  forces  or 
torques  at  a junction  equal  to  0,  i.e., 
f(t)  = 0 at  a junction 


(t)  = 0 at  a junction 

The  example  shown  in  Figure  4 illustrates  the  principles  discussed  for 
a translational  system.  This  example  represents  a translational  mechanical 
system  where  it  is  desired  to  solve  for  the  velocities  V2(t)  and  V^(t)  of 
the  system  when  an  external  force  is  applied  resulting  in  velocity  V^(t). 
Figure  5 represents  a mechanical  network  equivalent  of  Figure  4.  Equations 
(1)  and  (2)  result  when  the  velocity  measurements  .are  made  at  masses 
and  Mj. 


/•  r 

1 y (V''l)dt  + M2-5^+B(V2-V3)  + K2  j 


dt  = 0 


(1) 


r ‘^''3  r 

B (V3.V2)  + <2  j (V3-V2)  dt  . + K3  y^  V3dt 


(2) 
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Differentiating  each  side  with  respect  to  t results  in  the  set  of  second 
order  differential  equations  shown  by  Equations  (3)  and  (4). 


Ki  (V2'V^)  + + B (V2-V3)  + <2  (V2-V3)  = 0 (3) 

B (V3-V2)  + K2  (V3-V2)  + M3  V3  + K3  V3  = 0 (4) 


Figure  6 shows  an  electrical  network.  Equations  (5)  and  (6)  are  the 
differential  equations  generated  to  obtain  the  solution  for  voltages 
throughout  the  network. 
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(5) 

(6) 

and  are  identical  in  format  to  Equations  (3)  and  (4).  Figure  6 is 
the  electrical  analog  of  the  mechanical  network  shown  in  Figure  5. 

The  analog  relationships  generated  are  shown  in  Table  III. 
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Figure  6.  Electrical  Analogy  of  Mechanical  Example 
rr  ^2®2  ^ R'^®2  ■ ®3^  ^ I?  <®2  ' ®3^  " ° 

R-  ^®3  ■ ®2^  ^3®3  O'  ®3  " ° 


Table  III.  Mechanical/Electrical  Analogs  (Translational) 


Mechanical 

Electrical 

Relationship 

Force  (F) 

Current  (I) 

I = F 

Velocity  (V) 

Voltage  (V) 

V = V 

Spring  Stiffness  (K) 

Inductance  (L) 

L = 1/K 

Friction  (B) 

Resistance  (R) 

R = 1/B 

Mass  (M) 

Capacitance  (C) 

C = M 

The  dual  relationships  exist  for  tlic  rotatirnal  case  and  the  resulting 
conversion  is  shown  in  Table  IV. 


Table  IV.  Mechanical/Electrical  Analogs  (Rotational) 


Mechanical 

Electrical 

Relationship 

Torque  (T) 

Angular  Velocity  (W) 
Spring  Stiffness  (K) 
Friction  (B) 

Moment  of  Inertia  (J) 

Current  (I) 
Voltage  (V) 
Inductance  (L) 
Resistance  (R) 
Capacitance  (C) 

I = T 
V = W 
L = 1/K 
R = 1/B 
C = J 

The  analogous  quantities  in  the  table  are  used  directly  by  NELSIM  to 
obtain  an  electrical  analog  for  a given  mechanical  system.  The  non- 
electrical symbols  in  the  system  are  replaced  by  their  electrical  counter- 
parts and  thus  readied  for  analysis  using  a circuit  analysis  program.  The 
subroutines  utilized  to  achieve  this  are  discussed  in  Section  V. 
b.  Thermal/Electrical  Analogs 

Whenever  a temperature  gradient  exists  within  a system,  or  when  two 
systems  at  different  temperatures  are  brought  into  contact,  energy  is 
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transferred.  The  process  by  which  the  energy  transport  takes  place  is 
known  as  heat  transfer.  The  quantity  in  transit,  called  heat,  cannot  be 
measured  or  observed  directly,  but  the  effects  it  produces  are  amenable  to 
observation  and  measurement.  The  flow  of  heat,  as  in  the  performance  of 
work,  is  a process  by  which  the  internal  energy  of  a system  is  changed. 
Systems  which  are  primarily  concerned  with  the  transfer  of  heat  or  thermal 
energy  can  be  classified  as  thermal  systems. 

The  two  primary  elements  in  a thermal  system  are  the  portions  that 
act  as  heat  dissipators  (providing  thermal  resistance)  and  the  portions 
that  act  as  heat  reservoirs  (providing  thermal  capacitance).  These 
elements  are  connected  to  various  heat  sources  and  heat  sinks  (points  of 
varying  temperatures)  which  provide  the  potential  for  the  heat  flow  and 
are  related  to  the  rate  of  that  heat  flow  q(t)  as  shown  by  Equations 
(7)  and  (8). 

q(t)  = C,  ^ (7) 

q(t)  = the  rate  of  heat  flow  where  Ct  = C M,  Ct  = the  thermal 

T p T 

capacitance  of  the  mass,  Cp  = specific  heat,  M = mass  of  the 

body,  and  T = temperature  of  body. 

q(t)=l^T(t)  (8) 

where  Rj  is  the  thermal  resistance  and  may  be  expressed  as 

^ for  a region  of  lengtii  d with  uniform  cross-sectional 

area  A. 

As  seen  above,  the  form  of  the  equations  describing  the  response  of  a 
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thermal  system  depends  on  the  geometry,  configuration  and  type  of  the 
system  elements.  In  addition  to  these  requirements,  three  distinct  modes 
of  heat  transmission  are  recognized  and  must  be  specified.  The  three 
modes  are  conduction,  convection  and  radiation. 

The  basic  relation  for  heat  transfer  by  conduction  is  given  by 


dx  dt 


(9) 


where  = rate  of  heat  flow  by  conduction  in  a material, 

K = thermal  conductivity  of  the  material, 

A = area  of  the  section  through  which  heat  flows  by 
conduction,  to  be  measured  perpendicularly  to  the 
direction  of  heat  flow,  and 

^ = the  temperature  gradient  at  the  section. 


The  rate  of  heat  transfer  by  convectirn  between  a solid  surface  and 
a fluid  may  be  computed  by  the  following  relation: 


qc  = h^,AAT 


(10) 


where  q^  = rate  of  heat  transfer  by  convection, 

A = heat  transfer  area, 

AT  = difference  between  the  surface  temperature  T^  and 
a temperature  of  the  fluid  T„  at  some  specified 
location  (usually  far  away  from  the  surface),  and 
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= average  unit  thermal  conductance  (often  called  the 
surface  coefficient  of  heat  transfer  or  the  con- 
vection heat  transfer  coefficient). 

Using  the  above  equation,  the  thermal  conductance  for  convective  heat 
transfer  can  be  defined  as 

= TT^A  (11) 

and  the  thermal  resistance  to  convective  heat  transfer  R^,  which  is  equal 
to  the  reciprocal  of  the  conductance,  can  be  defined  as 

R = (12) 

The  rate  of  heat  transfer  by  radiation  between  two  bodies  is  given  by 
the  following  relation 

qr  = (13) 


where  q 

r 

o 


A 


1 


net  rate  of  heat  transfer  by  radiation, 

Stefan-Boltzmann  constant  (0.1714  x 10"®  Btu/hr 

ft^R^^). 

surface  area  of  body  1 , and 


^1_2  = 3 modulus  which  modifies  the  equation  for  perfect 

radiators  to  account  for  the  emissivities  and 
relative  geometries  of  the  actual  bodies. 

In  many  engineering  problems,  radiation  is  combined  with  other  modes 
of  heat  transfer.  The  solution  of  such  problems  can  often  by  simplified 
by  using  a thermal  conductance  K^,  or  a thermal  resistance  R^,  for  radi- 
ation. If  the  heat  transfer  by  radiation  is  written  as 


(14) 


the  conductance,  by  comparison  with  the  previous  equation,  is  given  by 


oA^fl-2(Ti-T2) 


(15) 


and  the  unit  thermal  conductance  for  radiation  h^  by 

\ = r^ 


(16) 


Ti-T^ 


where  is  any  convenient  reference  temperature. 

The  units  applicable  to  thermal  systems  that  will  be  accepted  by  the 
module  are  listed  in  Table  V. 


Table  V.  Thermal  Units 


Symbol 

Description 

Uni  t 

T 

Temperature 

op 

q 

Rate  of  Heat  Flow 

Btu  /sec 

«t 

Heat  Energy 

Btu 

t 

Time  in  Thermal  Circuit 

sec 

«t 

Thermal  Resistance 

(°F)(sec)(Btu) 

^t 

Thermal  Capacitance 

Btu/°F 

As  an  example  consider  Figure  7 which  is  used  to  illustrate  the  derivative 
of  a set  of  differential  equations  for  a thermal  system. 
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(a)  Themal  Satellite 
Problems 


(b)  Enlarged  Section  to  be  Studied 


Figure  7.  Thermal  System  Example 

In  the  thermal  system  example.  Figure  7 illustrates  an  orbiting  satellite 
with  one  outer  wall  facing  the  sun  and  one  facing  the  earth.  It  is  desired 
to  determine  the  temperature  within  its  four  adjacent  compartments  (the  two 
outside  temperatures  are  known).  If  the  compartment  walls  are  character- 
ized with  thermal  resistances  of  Rq,  , R2,  and  R^  and  each  compartment 
containing  thermal  capacitances  of  , C^.  and  and  temperatures  T-| , 
T2,  Tj  and  respectively,  the  thermal  analysis  model  is  given  by  the 
following  equations; 
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Analysis  of  the  electrical  network  shown  in  Figure  8 results  in  equation 
set  (19)  below  which  is  identical  in  form  to  Equation  set  (18). 


r 


^1  dT  = r (^-n-^l)  - ^ 

1 1 

1 1 

^4  dt^  " ^ ^®3'^'4^  ■ ^ ^®4"®out^ 


> 


(19) 


The  thermal  to  electrical  analogues  are  shown  in  Table  VI. 


Table  VI.  Thermal /Electrical  Analogues 


Thermal 

Electrical 

Relationship 

Temperature  (T) 

Heat  Flow  (Q) 

Thermal  Resistance  (R^) 
Thermal  Capacitance  (Cj) 

Voltage  (V) 

Current  (I) 

Electrical  Resistance  (R^) 
Electrical  Resistance  (C^) 

V = T 
I = Q 

Re  = Rt 
Ce  = Cj 
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The  practical  approach  to  thermal  problems  of  significant  complexity 
is  to  let  the  user  generate  the  thermal  analog  of  his  system  based  on 
required  physical  assumptions  and  approximations.  The  thermal  analog  can 
then  be  entered  into  the  NELSIM  module  in  terms  of  the  thermal  quantities 
given  in  Table  VI.  In  the  thermal  analogy  the  nodes  represent  thermal 
areas  connected  together  by  heat  paths  and  thermal  capacitances.  The 
thermal  analog  portion  of  NELSIM  accepts  this  type  of  input  and  converts 
it  to  electrical  analog  networks  . 

c.  Electro-Mechanical/Electrical  Analogs 

Electro-mechanical  systems  are  composed  of  both  electrical  and 
mechanical  devices.  The  interface  or  coupling  between  the  electrical  and 
mechanical  portions  of  the  system  is  effected  by  means  of  electro-mechanical 
transducers.  A transducer  is  a device  which  converts  energy  from  one  form 
to  another,  or  which  converts  a system  variable  from  one  form  to  another 
(e.g.,  a generator  converts  mechanical  energy  to  electrical  energy  and  a 
motor  converts  electrical  energy  to  mechanical  energy).  The  basis  for 
electro-mechanical  transducers  is  the  interaction  of  mechanical  forces  and 
magnetic  or  electrical  fields.  Thus,  the  unique  features  of  analysis  for 
electro-mechanical  systems  reduces  to  the  study  of  electro-mechanical 
coupling  through  magnetic  and/or  electric  fields. 

The  elements  of  an  electro-mechanical  system  can  be  visualized  as 
consisting  of  three  basic  types  as  illustrated  in  Figure  9:  electrical 

network  elements,  mechanical  system  elements,  and  electro-mechanical 
coupling  elements 


Electrical 

Electro- 

Mechanical 

Coupling 

Mechanical 

Network 

System 

Figure  9.  An  Electro  Mechanical  System 


The  differential  equations  describing  the  behavior  of  an  electro- 
mechanical system  are  based  on  the  assumption  of  the  validity  of  a lumped 
parameter  electro-mechanical  system.  In  the  context  used  here  lumped  para- 
meter systems  are  defined  as  follows:  the  electromagnetic  fields  are 
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quasistatic  and  electrical  terminal  properties  can  be  described  as  functions 
■|  of  a finite  number  of  electrical  variables.  Also,  the  mechanical  effects 

can  be  described  by  a finite  number  of  mechanical  variables.  Kirchoff's 
laws  can  then  be  written  for  the  electrical  parts  of  the  system  by  intro- 
ducing electro-mechanical  coupling  effects  through  the  terminal  relations 
of  the  coupling  system.  Similarly,  Newton's  second  law  and  continuity  of 
space  for  the  mechanical  parts  of  the  system  can  be  written,  including 
electro-mechanical  coupling  effects  in  the  terminal  relations  of  the 
coupling  system. 

The  NELSIM  program  allows  the  user  to  define  the  purely  electrical 
circuit  portion  of  the  system  in  terms  of  the  standard  electrical  elements 
(resistance,  capacitance  and  inductance),  the  purely  mechanical  portion  of 
the  system  in  terms  of  the  standard  mechanical  elements  (mass,  spring  and 
dampers),  and  the  coupling  terminals  in  terms  of  the  magnetic  flux  (or 
electrical  charge  for  an  electric  field  system). 

From  the  input  description  NELSIM  derives  the  electrical  equivalent  of 
the  electro-mechanical  system.  The  first  step  is  to  transform  the  mechan- 
ical portion  of  the  system  into  an  electrical  equivalent.  This  can  be  done 
utilizing  the  method  presented  in  Part  a.  of  this  section.  The  second  step 
is  to  connect  the  electrical  portion  of  the  system  to  the  generated  electrical 
analog  of  the  mechanical  portion  through  elements  and  equations  expressing 
the  coupling  between  the  two  portions.  To  illustrate,  consider  the 
galvanometer  sketch  provided  in  Figure  10.  A galvanometer  ideally  produces 
a deflection  that  is  proportional  to  an  electrical  input. 
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Figure  10.  Galvanometer  Schematic  Representation 

Figure  11  represents  the  equivalent  lumped  network  in  terms  of  its 
electrical  and  mechanical  parts 

ELECTRICAL  MECHANICAL 


Figure  11.  Lumped  Element  Galvanometer  Equivalent 

where  and  represent  the  moment  of  inertia  of  the  rotor  and  the  vis- 
cous damping  that  results  from  air  friction.  The  element  Kp  represents 
the  rotational  stiffness  produced  by  a spring  attached  to  the  rotor.  The 
constant  D is  defined  to  be  sta  where  6 is  the  magnetic  flux  density,  t 
is  the  total  length  of  the  conductor  and  a is  the  radius  of  the  rotor. 

The  two  dependent  sources  E^  and  provide  the  interface  between  the  two 
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systems  and  are  entered  directly  in  terms  of  the  elements  involved.  The 
electrical  analog  derived  by  NELSIM  is  shown  in  Figure  12.  As  can  be 
seen,  the  network  is  entirely  in  terms  of  electrical  quantities. 


Figure  12.  Galvanometer  Electrical  Analog 

Most  electro-mechanical  systems  can  be  modeled  in  terms  of  transfer 
functions  and  a user  has  the  option  of  entering  the  system  in  terms  of  a 
signal  flowgraph  (Section  III) . Models  with  general  usage,  however,  can 
be  easily  included  in  the  NELSIM  code.  To  illustrate  this  option  two 
models  have  been  built  into  NELSIM  and  are  discussed  below. 

(1)  Rate  Gyro  Model 

Rate  gyros  are  used  widely  in  any  application  requiring  tracking  and 
stabilization.  Aircraft  and  missiles  utilize  gyros  for  stabilization 
when  angular  deviations  are  noticed.  Rate  gyros  are  also  used  extensively 
in  conjunction  with  tracking  antennas  so  that  the  angular  velocity  of  the 
tracked  vehicle  can  be  measured.  A mathematical  block  diagram  of  a rate 
gyro  is  given  in  Figure  13. 
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Figure  13.  Block  Diagram  of  Rate  Gyro 

In  the  figure  represents  the  rate  input  to  the  gyro  and  e^  the  pro- 

portional output  voltage  generated.  The  K's  represent  amplifiers  within 
the  system  to  attain  the  desired  signal  levels  with  representing  the 
signal  generator  which  generates  an  electrical  signal  proportional  to  the 
displacement  angle  0^.  The  transfer  function  between  T2  and  T1  is  cal- 
culated as  follows: 

li  = K/C 
S(l+is) 

TTM- 

S(l+Js) 

when  simplified  the  above  equation  becomes 


T2  _ 

TT” 


1 + 


2 


The  above  represents  the  closed  loop  transfer  function  of  the  rate  gyro 
which  can  be  recognized  as  having  the  form 


F(S)  = 


1 + 


“n 


1 


n 
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where  5 is  the  damping  factor  and  the  natural  frequency.  Like  terms 
yield 

“n  ■ ^ ^ 

which  allows  the  user  to  adjust  both  of  these  quantities  by  adjusting  the 
various  gains. 

The  model  equivalent  for  NELSIM  input  is  shown  in  Figure  14  where 
PB  = 1/C  and  PC  = J.  The  feedback  gain  PF  has  been  added  for  further 
generalization  of  the  model  and  if  set  equal  to  one  will  result  in  the 
transfer  function  derived  above. 


Figure  14.  NELSIM  Equivalent  Rate  Qyro 

Solution  of  the  model  requires  breaking  the  equations  applicable  down 
to  a set  of  first  order  differential  equations  algebraically  coupled 
together.  Equation  (20)  gives  the  results  of  each  suimiing  junction. 
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PGl  = PIN-PA-P0UT-PF-PG2-PC 


(20) 

Solution  for  the  highest  order  derivatives  and  rearranging  (20) 
yields  Equation  set  (21). 

PGl  = PIN-PA-P0UT-PF-PG2-PC  \ 

PG2  = PGl-PB  ( (21) 

POUT  = PG2-PD  I 

A call  to  the  gyro  model  from  the  input  language  triggers  the  program  to 
generate  Equation  set  (21).  The  name  given  to  the  model  by  the  user  is 
added  to  each  variable  name.  For  example  if  the  model  is  referred  to  by 
the  user  as  T1 , parameter  PG7  becomes  PG7T1 . 

(2)  Accelerometer  Model 

The  accelerometer  is  the  basic  measuring  element  in  an  inertial 
navigation  and  guidance  system.  It  is  the  function  of  the  accelerometer 
to  provide  measurement  and  integration  of  linear  acceleration  to 
successfully  accomplish  the  functions  of  position  and  velocity  indication 
and  generation  of  proper  steering  signals.  Although  numerous  models  are 
available  the  double  integrating  accelerometer  shown  in  Figure  15  has 
been  implemented  into  NELSIM. 

In  the  model  a^  represents  the  input  acceleration,  e^  the  signal 
voltage  generated  and  0 the  angular  acceleration  of  the  motor  shaft 
resulting  from  the  application  of  e . The  double  integration  of  the 
angular  velocity  (o)  yields  the  angular  velocity  and  relative  angle 
(o)  of  the  shaft.  These  two  quantities  are  proportional  to  the  craft 


PG2  = PGl-PB'l 
POUT  = PG2-PD^ 
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velocity  (V)  and  distance  traveled  (X).  The  mathematical  formulation  is 
given  below. 


where: 


T3  _ 

S(1  + S S)  (1  + KgS)  + (K2K3gK^)S 


N 

1 + aS  + bS^  + N 


N = 


Inserting  T = K-ja^  and  applying  the  final  value  theorem  for  an  assured 
step  input  yields 

T3  = KK^a. 

in  the  steady  state  which  shows  the  angular  acceleration  of  the  shaft  to 
proportional  to  the  input  acceleration.  Another  relationship  exists  for 
T3,  namely 

T3  = Kg0S^ 

equating  the  two  yields 

KK^a^  = Kg0S^ 

To  determine  velocity  and  distance  changes  as  a function  of  a^  both  sides 
of  this  equation  are  integrated.  I 
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Digital  Pickoff 


Figure  15.  Accelerometer  Model 


The  model  progratimed  into  NELSIM  is  illustrated  in  Figure  16. 
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The  equations  obtained  at  the  sunning  junctions  are  shown  below. 


PA1  = PIN-PA-PA4 
PA3  = PA2-PD 


To  obtain  PA4  the  motor  transfer  function  block  between  PA3  and  PA4  is 
broken  down  into  a denominator  block  and  numerator  block  with  an 
additional  parameter  PA6  inserted  (Figure  17). 


Figure  17.  Breakdown  of  Motor  Transfer  Function  Block 


The  resulting  equations  for  the  motor  block  and  remainder  of  the  system 
are  shown  below. 


PA6  = PA3* 


1 

rrpF-s 


PA4  = PA6-PE-S 
PA5  . PA4^.^ 

POUT  = PA5-|- 

Multiplying  the  above  equations  out  and  solving  for  the  highest  order 
derivatives  yield  Equation  set  (22). 


* 

^ 
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PAl  = PIN-PA-PA4 

PA2  = (PA1-PB-PA2)/PC 

PA3  = PA2-PD 

f’A6  = (PA3-PA6)/PF 

PA4  = PA6-PE 

PA5  = PA4/PG 

POUT  = PA5 


(22) 


Rearranging  Equation  set  (22)  and  breaking  down  n nurnber  equations  to  n 
first  order  differential  equations  yields  the  set  of  Equations  (23)  shown 
below. 


PAl  - PIN-PA-PA4 
PA2  = PA7 

PA7  = (PA1.PB-PA7)/PC 

PA3  = PA2-PD  ^ (23) 

PA6  = (PA3-PA6)/PF 
PA4  = PE-(PA3-PA6)/PF 
PA5  = PA4/PG 
POUT  = PAS 

A call  to  the  accelerometer  model  triggers  the  program  to  generate 
Equation  set  (23).  The  program  adds  to  each  parameter  involved  the  name 
assigned  to  the  model  by  the  user.  For  example,  referral  to  the  model  as 
Ml  will  result  in  the  letters  Ml  being  added  to  each  parameter  (i.e., 
PA4M1 ) . 
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d.  tlectro-Optical/Electric  Analogs 

The  majority  of  photodectors  in  use  today  are  photo  emissive  diodes, 
photo  conductive  diodes,  or  PIN  and  PN  photo  diodes.  Photo  emissive  diodes 
typically  have  low  quantum  efficiencies,  while  photo  conductive  diodes  have 
low  cut-off  frequencies.  These  properties  have  tended  to  push  the  using 
community  towards  photo  diodes  except  where  special  applications  make  the 
other  devices  specifically  attractive.  A direct  result  of  this  trend  has 
produced  a large  quantity  of  photo  diode  data,  devices  and  literature,  and 
lesser  quantities  of  the  other  device  types. 

The  NELSIM  translator  has  been  constructed  to  handle  photo  diode  type 
devices,  simulating  the  AC  characteristics  of  a device  with  an  equivalent 
electric  circuit.  Circuit  representation  is  shown  in  Figure  18. 


Figure  18.  Eqvivalent  Circuit  of  a Photodiode 

Referring  to  the  figure,  Ip^  represents  the  AC  photocurrent  generated 
by  the  photo  diode  when  excited  by  a light  energy  source.  , R^  , and 
C^  are  the  equivalent  resistance  and  capacitance  that  characterize  the 
device.  Table  VII  presents  performance  characteristics  for  widely 
used  photo  diodes. 

The  photocurrent  can  be  determined  from  light  energy  time  histories 
using  perfect- square  law  detector  formulation.  Typically,  the  photo 
current  is  expressed  as  some  function  of  the  quantum  efficiency  and  average 
nurrtier  of  incident  photons  per  unit. 
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time,  such  as 


IpH  ' (nqP/hv) 


where , 


ri  = quantum  efficiency 

P/hv  = average  number  of  incident  photons 
per  unit  time 

Ipi^/q  = average  number  fps  unit  time  of 

electrons  emitted  from  the  photocuthode 


Since  the  device  response  is  dependent  upon  the  frequency  of  the  energy 
source  it  is  best  to  construct  a circuit  for  frequency  bands  and  sum  the 
produced  current  linearly,  or  to  restrict  the  analysis  to  a certain 
frequency  band.  Figure  19  details  the  effect  of  frequency  on  quantum 
efficiency  for  various  photodiode  devices. 


Figure  19.  Wavelength  Dependence  of  Quantum  Efficiency  and 

Responsitivlty  for  Several  High  Speed  Photodiodes. 
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j e.  Units  and  Scaling 

The  following  paragraphs  document  steps  taken  to  implement  units  and 
a scaling  algorithm  into  NELSIM.  In  order  to  allow  NELSIM  to  be  completely 
compatible  with  SCEPTRE  in  the  flexibility  of  its  input  language  the 
scaling  module  is  semi  automated  and  requires  user  interaction. 

(1)  Units 

SCEPTRE  allows  any  set  of  parameter  units  as  long  as  the  element 
magnitudes  with  respect  to  these  are  consistent.  To  keep  the  two 
input  languages  consistent,  the  same  is  true  of  NELSIM.  The  user  is 
allowed  to  enter  a problem  in  any  set  of  units  he  chooses  and  is  responsible 
for  compatibility  throughout  his  system. 

(2)  Scaling  Guidelines  and  Automation 

One  of  the  problems  encountered  in  system  analysis  utilizing  numerical 
r 

integration  is  that  of  maintaining  a consistent  set  of  parameter  units  such 
that  the  magnitudes  of  all  the  solution  variables  are  consistent  with  the 
specified  error  criteria.  Numerical  inaccuracies  result  and  excessive 
solution  time  is  spent  when  units  are  not  carefully  chosen  due  to  the 
difficulties  in  meeting  specified  error  criteria.  Previous  experience 
in  the  system  analysis  area  and  application  of  SCEPTRE  and  other  codes  to 
large  scale  circuit  and  system  problems  has  shown  that  it  is  desirable  to 

_3 

maintain  the  magnitude  of  the  solution  variables  in  the  range  between  10 

3 

and  10  . For  example,  voltage  and  current  are  the  two  solution  variables 
used  by  SCEPTRE  and  it  is  desirable  to  keep  these  variables  within  the 
limits  shown  below. 

10"^  < |Vl  ^ llO^ 

10"^  i |i|  1 
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The  desired  magnitude  ranges  can  be  accomplished  through  knowledge  of  the 
basic  laws  governing  a particular  system  and  a knowledge  of  the  ranges 
expected.  If  the  expected  units  of  the  solution  variables  and  time  are 
known,  the  system  element  units  can  be  scaled  to  yield  those  units.  The 
above  is  exemplified  through  a high-speed  transistorized  circuit  in  which 
the  response  is  known  to  be  in  terms  of  volts,  mi  Hi  amps  and  nanoseconds. 
The  equations  below  illustrate  the  applicable  relationships  between  state 
variables  and  electrical  elements. 


dt 

^ ^ "37 


Where:  V = Voltage 
I = Current 
R = Resistance 
C = Capacitance 
L = I nductance 


To  obtain  the  response  in  terms  of  volts,  milliamps  and  nanoseconds  the 
elements  are  scaled  as  follows: 


n _ y.  = = -i--y  = 10^  (kilohms) 

I ma  iq-3 

r - T - ma  ^ = 10'^^  (Picofarads) 

C - I ^ - ma.  T 


L = V 


^ = Volts  . 


nsec 

ma 


10”^  (microhenries) 


Definition  of  the  network  elements  in  terms  of  these  units  brings  the 
solution  variables  closer  to  their  desired  range. 
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The  same  procedure  is  applicable  to  other  systems.  The  equations 
shown  below  present  the  relationships  necessary  for  mechanical  and  thermal 
systems. 


MECHANICAL  SYSTEMS 


M = 


- i ^ 

■ V dt 


Where:  V = Velocity 

F = Force 
M = Mass 

K = Spring  Constant 

B = Viscocity/ 
Friction 


THERMAL  SYSTEM 

- _ dt  Where:  Q = rate  of  heat 

^ dT  flow 

y C = Thermal  Capacitance 

= Q 

^ R = Thermal  resistance 

T = Temperature 

(3)  Automated  Scaling 

An  algorithm  was  implemented  within  the  NELSIM  program  to  achieve 
scaling  of  a network.  The  program  generates  an  electrical  analog  from 
various  types  of  system  inputs.  The  scaling  algorithm  scales  all  the 
elements  to  magnitudes  consistent  with  user  provided  guidelines.  The 
option  requires  the  user  to  input  the  expected  units  of  the  solution 
variable  and  time.  Table  VIII  shows  the  analogies  between  the  various 
systems  involved.  Both  mechanical  and  thermal  quantities  are  translated 
into  the  electrical  quantities  shown  in  the  table,  such  that  in  the  trans- 
lated network  the  only  solution  variables  are  electrical.  The  user  will 
enter  a problem  in  terms  of  the  standard  units  familiar  to  the  field 
involved.  Table  IX  illustrates  the  units  that  are  expected  if  the  scale 
option  is  utilized. 
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TABLE  VIII.  SYSTEM  ANALOGIES 


SYSTEM 

FUNCTION 

ELECTRICAL 

MECHANICAL 

THERMAL 

Translational 

Rotational 

SOLUTION 

VARIABLES 

Voltage 

Velocity 

Angular  Velocity 

Temperature 

Current 

Force 

Torque 

Heat  Flow 

TIME 

Time 

Time 

Time 

Time 

ELEMENTS 

Resistance 

Friction 

Friction 

Thermal 

Resistance 

Capacitance 

Mass 

Moment  of  Inertia 

Thermal 

Capacitance 

Inductance 

Spring  Stiffness 

Spring  Stiffness 

TABLE  IX.  STANDARD  SYSTEM  UNITS 


PHYSICAL  MfniUM 


Electrical 


ELEMENTS 


Resistance 

Capacita^-p 

Inductance 


Mass 

Spring 

Stiffness 

Friction 

Moment  of 
Inertia 


Angular 

Velocity 


SYMBOL  I STANDARD 


Ohms 


UNITS 

ENGLISH  BTU 


Kilogram 

Slug 

Newton/ 

meter 

lb/ ft 

n-sec/m 

Ib/ft 

nm/  , 

(rad/sec^) 

ft  Tb/  , 

(rad/sec^) 

rad/sec 

rads/sec 

TABLE  IX.  (CONTINUED) 


PHYSICAL  MEDIUM 

ELEMENTS 

SYMBOL 

STANDARD 

MKS 

ENGLISH 

BTU 

Thermal 

Temperature 

T 

°F 

Heat  flow 
rate 

Q 

Btu  /sec 

Capaci tance 

c 

Btu /°F 

Reiistance 

R 

8tu/°F/ 

lb 

The  user  must  indicate  the  scaling  to  be  performed  by  entering  the  expected 
magnitudes  of  the  solution  variables  and  time.  For  example  consider  the 
electronic  network  discussed  earlier  entered  in  terms  of  the  standard 
elements  shown  in  Table  VIII.  To  obtain  the  solution  in  terms  of  volts, 
mil li amps  and  nanoseconds  the  network  must  be  scaled.  The  user  input 
would  consist  of  the  card  shown  below  which  states  that  the  scaling  option 
is  desired 

SCALE  OPTION,  V0LTAGE=1 , CURRENT=10E-3,  TIME=10E-9 
such  that  the  voltage  is  to  remain  in  terms  of  volts  but  current  and  time 
are  to  be  scaled  to  milliamps  and  nanoseconds  , respectively.  Use  of  this 
scaling  option  for  all  system  components  will  maintain  a consistent  set  of 
variable  values. 

The  output  of  the  program  consists  of  two  listings.  The  first  is  a 
listing  of  the  analog  network  generated  with  no  scaling  applied.  The 
second  consists  of  the  analog  network  with  all  constant  elements  scaled  and 
a list  of  the  scale  factors  used  on  the  elements.  The  user  is  then 
required  to  scale  all  non-constant  element  values  to  complete  the  scaling 
task.  This  requirement  is  a result  of  the  flexibility  of  the  input 
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language  given  tn  the  user  to  keep  the  program  consistent  with  SCEPTRE. 
Due  to  this  flexibility  the  user  has  the  ability  to  camouflage  non- 
uni tl  ess  quantities,  as  will  be  illustrated.  One  of  the  allowed  input 
formats  of  NELSIM  is  the  definition  of  an  element  as  a Fortran 
function  subroutine.  The  subroutine  may  contain  non-unitless  quantities 
which  would  have  to  be  scaled.  Yet  the  program  cannot  possibly  determine 
the  units  of  these  parameters  and  any  attempt  at  scaling  these  could 
result  in  gross  error.  Consider  the  case  of  a voltage  dependent  resistor 
entered  as  follows: 

R1,1-2=FUNCTI0N  R(ET) 

where  the  Function  R is  described  as 

FUNCTION  R(X) 

R = 10. *X 

RETURN 

END 

The  constant  10.  is  in  units  of  — and  any  scaling  of  current  would 
demand  a scaling  of  this  constant.  The  program,  however,  has  no  access  to 
the  function  and  no  way  oF  knowing  about  the  constant.  For  this  reason, 
the  task  of  altering  the  function  will  be  left  to  the  user. 

A completely  automatic  scaling  algorithm  would  demand  a severe 
reduction  in  the  flexibility  of  the  input  language  which  would  in  turn 
defeat  the  purpose  of  keeping  the  program  compatible  with  transient 
analysis  programs. 
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SECTION  III 


DIFFERENTIAL  EQUATIONS 

Section  II  points  out  that  differential  equations  are  a very  general 
form  of  system  definition.  This  section  presents  differential  equation 
capabilities  which  are  built  into  NELSIM. 

The  NonElectriv-al  Languages  Simulation  Module  accepts  coupled  sets 
of  nonlinear  (or  linear)  differential  equations  of  the  form  given  by 
equation  (24). 
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Equation  (24)  is  a general  set  of  m nonlinear  differential  equations  in 
m unknowns,  through  and  with  f^rcing  functions  through  F^.  For 
convenience,  the  equations  are  represented  in  matrix  form,  but  are  input 

d"x,’ 

input  into  NELSIM  individually.  The  term  — — is  the  nth  derivative 

dt" 

of  the  variable  raised  to  the  ith  power  where  i can  vary  from  1 to  i 
and  is  not  necessarily  the  same  for  all  derivative  orders.  For  the 
systems  considered  in  Section  II,  the  highest  order  of  derivative  n 
generally  encountered  is  two.  The  coefficients  a-ji,  b^-j,  etc.  may  be 
time  varying  or  nonlinear.  It  will  be  possible  to  enter  the  coefficients 
into  NELSIM  as  constants,  SCEPTRE  type  defined  parameters,  or  algebraic 
expressions,  e.g.,  the  form  (3t  - t + 1). 

NELSIM  automatically  takes  equations  of  the  form  of  Equation  (24) 
and  separate  them  into  first  order  differential  equations.  It  should  be 
realized  that  to  determine  the  transient  response  of  the  system  represented 
by  Equation  (24),  it  will  be  necessary  to  integrate  a number  of  equations 
equal  to  the  total  number  of  derivatives  up  to  n derivatives  for  each 
unknown  for  a maximum  of  m times  n equations  to  integrate.  The  maximum 
number  of  first  order  differential  equations  currently  allowed  in  SCEPTRE 
for  example  is  100,  and  NELSIM  warns  the  user  if  over  100  equations  are 
generated.  If  the  SCEPTRE  limit  were  increased,  NELSIM  could  be  easily 
modified  to  acconmodate  the  increase.  Each  nth  order  equation  must  be 
translated  into  n first  order  differential  equations  as  shown  below. 
Equation  (24)  is  broken  down  to 
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A ./^  AS 

dt"  V « • dt  • 


dx 

dt’ 


which  is  rewritten  as  Equation  set  (25). 


dx  , 
m- 1 _ „ 

"3t~- 


Vl 


..  C) 


\ 

> 

, 


(25) 


As  an  example  consider  the  two  second  order  differential  Equations  26  and 

27  derived  for  the  mechanical  system  of  section  I la  and  repeated  here  for 
convenience. 


Ki(V2-Vi)  + M2V2  + B(V2-V3)  + K2(V2-V3)  = 0 (26) 

B(V3-V2)  . K2  (V3-V2)  . M3V3  . K3V3  = 0 (27J 

If  entered  into  NELSIM  according  to  the  proper  format,  the  program  will 

proceed  to  solve  each  equation  for  the  highest  order  derivative  and  from 
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these  generate  a set  of  four  first  order  differential  equations  as  shown 
below. 

Solving  for  highest  order  derivatives  yields: 


V2  =[Ki(V''i)  ■ ^^V'^3^  “ ’^2^V''3^]  /^2 

(28) 

V3  =[Kj(V3-V2)  - 

(29) 

Breaking  (28)  and  (29)  into  four  first  order  differential 
yields: 

equations 

^2  ^^2A 

(30) 

^2A  = ■''2  ■ - ®<''2-*3>  - ^2<''2-''3>]  '"2 

(31) 

^3  " '^3A 

(32) 

^3A'"  '^*3  " ['*'2^'^3‘'^2^  ■ ®^^3'^2^  " ^3  ] "^^3 

(33) 

As  shown,  two  new  variables  are  generated  in  the  process  {^2^ 
Substitution  of  these  two  variables  into  Equations  (31)  and  (32)  yields 
the  set  of  equations  shown  below. 

^2A  " ^2^ 

(34) 

"[’^l^V'^l^  ■ ®^''2A"''3A^  ■ *^2^''2'''3^]  ^^2 

(35) 

^3  " '^3A 

(36) 

''3A'[MVV-8(W-  V3]'"3 

(37) 
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SECTION  IV 
TRANSFER  FUNCTIONS 

This  section  presents  the  transfer  function  capabilities  which  are 
built  into  NELSIM.  A system  defined  in  terms  of  linear,  time  invariant, 
ordinary  differential  equations  can  be  represented  by  a transfer  function 
which  relates  the  system  output  to  input  in  the  Laplace  transform  domain 
using  rational  polynomials  in  the  Laplace  transform  variable  S.  Converse- 
ly, given  the  transfer  function  of  a real  system  as  a rational  polynomial 
in  S,  it  is  always  possible  to  derive  a set  of  first  order  differential 
equations  with  time  as  the  independent  variable  whose  solution  is  the 
transient  response  of  the  system  for  a given  input.  The  purpose  of  the 
transfer  function  section  of  NELSIM  is  to  provide  the  capability  to  input 
a connected  set  of  block  diagrams  in  which  the  individual  blocks  are 
rational  polynomials  in  S.  The  program  automatically  derives  the  appro- 
priate first  order  differential  equations  . The  transient  analysis 
program  can  then  integrate  the  differential  equations  to  provide  the  desired 
transient  output. 

For  a transfer  function  of  the  form  given  by  Equation  (38) 


(38) 


where  C(s)  is  the  Laplace  transform  of  the  output  C(t)  and  R(s)  is  the 
Laplace  transform  of  the  forcing  function  R(t),  the  following  state 
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variable  differential  Equations  (39)  and  single  algebraic  equation  (40) 
define  the  transient  response. 


dt  “ ^2 

dx2 

W 


dx  , 
n-1  _ 

dt  ■ ^n 


L 

dt  b. 


R(t)  - b.x^-i+1 


i = l 


C(t)  ■ ^ 
i=0 


^Vi+1 


/ 


(39) 


(40) 


The  set  of  Equations  (39)  and  (40)  are  developed  automatically  by  NELSIM 
for  each  transfer  function  entered.  The  interface  between  transfer 
functions  is  taken  care  of  by  the  nodal  connections. 

As  an  example  of  the  transfer  function  capability  in  NELSIM,  consider 
the  unity  feedback  system  shown  in  Figure  20. 
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Figure  20.  Unity  Feedback  System  Example 

From  the  input,  the  program  derives  equations  at  each  node  as  follows: 
P2  = PI  - POUT 
P3  = (P2  - 1.2  P3  - P3)/.2 
POUT  = 3P3  + P3 

These  equations  are  then  broken  down  to  first  order  aquations  if  necessary 
as  shown  in  Equation  set  (41)  and  the  proper  variable  substitutions  made. 
The  generated  set  of  equations  acceptable  to  SCEPTRE  are  listed  below. 

P2  = PI  - POUT  \ 

P3  = P3A  / 

P3A  = P3  = P3B  / ("*1) 

P3B  = P3  = (P2  - 1.2  P3B  - P3A)/.2  1 

POUT  = 3 P3A  + P3  / 

The  extra  parameters  generated  by  NELSIM  (P3A  and  P3B)  are  unique  variables 
whose  names  are  generated  from  the  original  parameter  name  (i.e.  , P3). 
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SECTION  V 


I 

PROGRAM  CONFIGURATION 

The  NonElectri cal  Languages  Simulation  Module  consists  of  a 
FORTRAN  program  which  will  accept  as  input  descriptions  of  systems  other 
than  electrical  and  converts  these  to  a form  acceptable  to  the  SCEPTRE 
analysis  program.  The  output  of  this  module  is  in  the  form  of  a printed 
listing  and/or  deck  of  punched  cards  ready  for  immediate  analysis  by 
SCEPTRE.  Figure  21  below  illustrates  the  main  function  of  the  module. 


SYSTEM 

DESCRIPTION 


NELSIM 


INPUT  DECK  READY 
FOR  SCEPTRE  INPUT 


Figure  21.  Main  Function  of  Non-Electrical 
Languages  Simulation  Module 

The  input  language  of  NELSIM  is  format  free  and  user  oriented  with  a 
syntax  similar  to  that  of  SCEPTRE.  The  input  types  allowed  for  NELSIM 
are  illustrated  in  Figure  2 and  discussed  in  Sections  II  through  IV.  This 
section  documents  the  configuration  of  the  program  and  subroutines 
utilized  to  perform  the  functions  described  in  the  past  sections. 

The  program  is  made  up  of  three  modules,  an  input  processor, 
translator  and  output  processor.  The  input  processor  reads  in  the  input 
data,  determines  the  type  of  translation  needed,  and  stores  and  prepares 
the  information  for  the  translator.  The  translator  performs  all  necessary 
transformation  calculations  and  scaling.  It  then  stores  the  information 
in  a manner  easily  accesible  to  the  output  processor. 
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The  output  processor  accesses  the  information  stored  by  the  system 
translator  and  outputs  it  in  a format  compatible  with  the  SCEPTRE  input 
languages.  The  information  can  be  punched  on  cards  at  the  option  of  the 
user.  The  input  is  allowed  in  terms  of  combinations  of  blocks  containing 
either  differential  equations,  transfer  functions  or  functional  elements 
of  the  various  disciplines  coupled  by  nodal  connections  or  algebraic 
equations.  The  program  translator  treats  each  block  individually  per- 
forming the  appropriate  conversion  and  then  connects  them  as  defined  by 
the  user.  A flow  chart  of  the  main  functions  of  the  program  is  provided 
in  Figure  22. 


Figure  22.  NELSIM  Major  Functions 

The  translated  output  from  NELSIM  is  configured  to  produce  information 
of  a form  that  can  be  used  directly  as  input  into  SCEPTRE.  In  order  to 
allow  alteration  of  the  output  format  for  use  with  other  programs,  the 
routines  that  require  change  have  been  placed  in  the  output  processor  module 
and  have  been  given  subprogram  names  that  end  in  "OUT".  These  routines  are: 
ELOUT,  DPOUT , RCOUT,  ICOUT,  and  MODOUT.  A description  of  these  routines 
is  contained  in  Table  X. 
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Only  the  main  subroutines  and  major  flow  of  the  program  have  been 


I 


illustrated.  Table  X below  contains  a list  of  these  subroutines  and 
their  particular  functions.  A complete  list  of  all  subroutines  utilized 
by  NELSIM  is  provided  in  Table  X. 


Table  X.  Major  NELSIM  Routines  and  Functions 


ROUTINE  NAME 

FUNCTIONS 

DRIVER 

Main  program,  controls  entire  execution  of  program  and 
contains  calls  to  all  modules 

GTCRD 

Get  a card  and  determine  heading  or  subheading  type 

ELMSET 

Determine  starting  locations  of  elements,  defined 

DPSET 

parameters,  models,  output  and  set  initial  conditions 

MODSET 

OUTSET 

ICSET 

of  variables  for  each  system  entered 

ELMPRS 

Process  and  store  element  name,  nodal  connection  and 
value  codes 

DPPRS 

Defined  parameter  processor  to  store  defined  parameters, 
their  locations  and  value  codes 

OUTPRS 

Output  processor  to  process  and  store  all  outputs 

ICPRS 

Initial  condition  processor 

RCPRS 

Run  control  processor,  interprets  and  stores  all  run 
controls  information 

FNPRS 

Function  processor,  processes  and  outputs  all  tabual  or 
data  as  well  as  equation  data 

DET 

Driver  for  differential  equation  translator 

TFT 

Driver  for  transfer  function  translator 

SNODES 

Reorder  and  rename  nodal  connections  as  necessary 

CHANGE 

Translate  mechanical  elements  to  electrical  analogs 

mmni 

Translate  thermal  elements  to  electrical  analogs 

Output  all  elements  under  the  'Elements'  subheading 

ELWRITE 

Write  out  elements  and  values  onto  temporary  tape  prior 
to  final  output  by  DRIVER. 
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Table  X.  Major  NELSIM  Routines  and  Functions  - Continued 


ROUTINE  NAME 

FUNCTIONS 

EQFCHG 

Change  contents  of  equation,  function  or  expression  from 
user  provided  quantities  to  electrical  analog  quantities 

DPOUT 

Output  defined  parameters  under  'DEFINED  PARAMETERS' 
subheadi ng 

DPWRITE 

Write  out  defined  parameters  onto  temporary  tape  prior 
to  final  output  by  DRIVER 

RCOUT 

Output  all  run  control  information 

I CO  LIT 

Output  all  initial  conditions 

OUTPTO 

Output  all  outputs  of  system 

OUTWRIT 

Write  outputs  on  temporary  tape  prior  to  final  output 
by  DRIVER 

SCALE 

Scale  element  values  to  units  provided  by  user 

MODFUNC 

Write  out  functions  necessary  for  models  used 

MODOUT 

Driver  to  output  models  desired  in  defined  parameter 
format 

GYRO 

Subroutine  containing  first  order  differential  equations 
describing  the  gyro  built  in  model 

MODWRIT 

Write  model  information  on  temporary  tape  prior  to  actual 
output  by  DRIVER 

ACCEL 

subroutines  containing  first  order  differential  equations 
necessary  to  describe  the  built  in  accelerometer  model 

Each  nxjdule  is  discussed  in  detail  in  the  next  paragraphs. 


Figure  23  illustrates  the  flow  of  the  input  processor.  The  functional 
subheading  processor  is  shown  in  Figure  24  and  stores  all  the  information 
regarding  the  functional  subheading,  for  example  ELMPRS  stores  all  ele- 
ment related  information  such  as  element  type,  nodal  connections  and 
element  value. 

The  major  subroutines  utilized  by  the  analog  translator  are  shown  in 
Figure  25.  These  subroutines  generate  electrical  analogs  from  input 
thermal,  mechanical,  electro-mechanical  and  electro-optical  system 
descriptions  as  defined  in  Section  II. 


51 


CALI  lAMSFIR  ruNCTlOM 
TMUSlATOH 


«A0  l«klN 
HtAOtHG  TVn 


« I UUTPUT 

i I PROUSSOd  I 


TMUS^ER 

ruNCTions. 


VES  ^iFrFBfwTIAP* 
N.  EOUATlONSv 


Ml  »«CMANlCAl 
EiAfiS 


sn  THEflmL  a*Gsl 


CM.L  DlPrtPCNTIAL 
COuATIONS 
translator 


YES  ^ ELtaRO-'' 
\ »€04ANICAL. 


SET  t/m  R.AGS 


^ ANALOG  ^ 

.:«SCRIPTION. 


YES  ^ ELEaRO- 

\ optical 


sn  e/o  aAGs 


PiAiaiOPlAL 

.SUBHEAOINC. 


CAU  SUVCAOING  PROCESSOR 


Figure  23.  NELSIM  Input  Processor 
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Figure  24.  Function  Subheading  Processor 
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Figure  25.  Translator  Main  Routines 
5 


When  differential  equations  are  entered  the  program  calls  on  the 


Differential  Equation  Translator  (DET)  which  performs  the  functions 
documented  in  Section  III.  The  functional  flow  of  the  DET  module  is  as 
follows : 

1.  Locate  all  delimiters  within  equation. 

2.  Break  equation  into  major  parts. 

3.  Isolate  parts  containing  derivative  terms  on  left  side  of 
equal  sign  and  algebraically  move  remaining  terms  to  right 
of  equal  sign. 

4.  Isolate  highest  derivative  term  and  algebraically  move 
remaining  derivative  terms  to  right  of  equal  sign. 

5.  Isolate  highest  derivative  term  from  its  coefficient, 
moving  coefficient  algebraically  to  right  side  of  the 
equation. 

6.  Examine  highest  derivative  order  and  convert  order 
equation  to  N 1^^  order  equations. 

Figure  26  depicts  the  flow  for  a sample  equation. 

The  Transfer  Function  Translator  (TFT)  is  utilized  when  a set  of 
transfer  functions  are  input  into  the  program.  The  TFT  generates  the  set 
of  first  order  differential  equations  depicting  the  system.  The  flow  of 
the  TFT  module  is  illustrated  in  Figure  27. 

The  output  processor  is  illustrated  in  Figure  28.  The  functions 
such  as  tables  and  equations  are  output  by  the  program  first,  then  all 
the  elements  are  output  followed  by  all  the  defined  parameters.  The 
requested  outputs  are  then  written  out.  Output  of  initial  conditions, 
run  control  information  and  rerun  description  finishes  the  output  function. 
The  two  output  options  are  then  executed,  namely  if  the  punch  option  is 
requested,  the  deck  is  punched  and  if  the  scale  option  is  requested,  the 
scaled  output  is  listed  out. 
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DIFFERENTIAL  EQUATION  TRANSLATOR 


I 

1.  GIVEN  THE  FOLLOWING  EQUATION 
B*(D1V3-D1V2)  +K2*(V3-V2)+M3*D2V3+K3*V3=F1*TI 

2.  COMPRESSING  OUT  ALL  BLANKS  YIELDS 

+B*( D1 V3-D1 V2 )+K2*( V3-V2)+M3*D2V3+K3*V3=F1 *T1 

3.  BREAKING  THE  EQUATION  INTO  PARTS  YIELDS 

+B*(D1V3-D1V2) 

+K2*(V3-V2) 

+M3*D2V3 

+K3*V3 

=F1*T1 

4.  REARRANGING  THE  PARTS  YIELDS 

+B*(D1V3-D1V2) 

+M3*D2V3 

=F1*TI  -K2*(V3-V2)  -K3*V3 

5.  ISOLATING  ONLY  THE  HIGHEST  DERIVATIVE  YIELDS 
+M3*D2V3 

=FI*T1  -K2*(V3-V2)  -K3*V3  -B*(DI V3-D1V2) 

6.  SEPARATING  THE  DERIVATIVE  TERM  AND 
COMPRESSING  OUT  THE  BLANKS  YIELDS 

D2V3=(F1*T1-K2*(V3-V2)-K3*V3-B*(D1V3-D1V2))/M3 

7.  TRANSFORMING  THE  EQUATION  YIELDS 

DPV3A=(PF1*PT1-PK2*(PV3-PV2)-PK3*PV3-PB*(PV3A-PV2A))/PM3 

DPV3=PV3A 


Figure  26.  DET  Functional  Flow 


J 


56 


START 


Figure  27.  Transfer  Function  Translator  Functional  Flow 
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Figure  28.  Output  Processor  Routines 
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